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A UNIFIED FORM OF LAMBERT'S THEOREM 

E. R. Lancaster 
R. G. Blanchard 


ABSTRACT 

A unified form of Lambert's theorem is presented which is valid for ellip- 
tic, hyperbolic, and parabolic orbits. The key idea involves the selection of an 
independent variable x and a parameter q such that the normalized time of flight 
T is a single-valued function of x for each value of q. The parameter q depends 
only upon known quantities. If T is less than the time of flight for one complete 
revolution, it is a monotonic function of x for each q, making possible the con- 
struction of a simple algorithm for finding x, given T and q. Detailed sketches 
are given for T(x, q) and formulas developed for the velocity vectors at the in- 
itial and final times. Also included is a careful derivation of the classical form 
of Lambert's equations, including the multi-revolution cases. 
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A UNIFIED FORM OF LAMBERT'S THEOREM 


INTRODUCTION 

Lambert's problem, as it arises in most applications, is concerned with the 
determination of an orbit from two position vectors and the time of flight. It has 
important applications in the areas of rendezvous, targeting, guidance, and 
preliminary orbit determination. In this paper a unified form of Lambert's 
theorem will be presented which is valid for elliptic, hyperbolic, and parabolic 
orbits. 

The key idea involves the selection of an independent variable x and a param- 
eter q such that the normalized time of flight is a single valued function of x for 
each value of q. The parameter q depends only upon known quantities. The 
problem then is to find x for given values of q and the time of flight. If the time 
of flight T is less than that for one complete revolution, T is a monotonic func- 
tion of x for each value of q. Thus it is an easy task to design an algorithm for 
finding x. For multirevolution cases T(x) has a single well-defined minimum 
for each q„ 

This idea was presented in a previous paper [1], where a unified formula was 
given for the computation of T from x and q. Detailed sketches were given for 
T(x, q) and a simple formula developed for the time derivative of the magnitude 
of the radius vector at the initial time in terms of x and given quantities. 

The present paper expands upon the previous one by 

1. giving complete derivations which were only sketched before, 

2. giving a careful derivation of the classical form of Lambert's equations, 
including the multirevolution cases, 

3. deriving a number of useful auxiliary formulas, e.g., for the semilatus 
rectum and for the velocity vectors at the initial and final times in terms 
of the two given position vectors. 


THE PROBLEM 

Suppose a particle in a gravitational inverse-square central force field has 
distances r x and r 2 from the center of attraction at times t x and t 2 . Let c be 
the distance and 6 the central angle between the positions of the particle at the 
two times, where 0 <6<2tt. 
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Lambert's problem is that of finding the semimajor axis or some related 
quantity for the orbit of the particle, given t x , r ls t 2 , r 2 , and 9. Having solved 
Lambert's problem, other quantities associated with the orbit are easily found, 
as will be later discussed. Using the law of cosines, we can express c in terms 
of r x , r 2 , and 6: 


c 2 = r x 2 + r 2 2 “ 2r i r 2 cos 9 


( 1 ) 


We define 

G - universal gravitational constant 
M = mass of attracting body 

/x = GM 

a = semimajor axis of transfer orbit 
e = eccentricity of transfer orbit. 

We will follow the common sign convention for a, i.e., a > 0 for elliptic 
orbits and a < 0 for hyperbolic orbits. Definitions of other symbols will be 
given as they are introduced. 


THE CLASSICAL FORM OF LAMBERT'S EQUATIONS 

The path of a particle in an inverse-square central force field is an ellipse, 
parabola, or hyperbola. With origin at the center of attraction, we have, for 
elliptic motion 


(l - e cos , 

(2) 

i (l - e cos 0 2 ) , 

(3) 

= (pj - e sin cp 1 , 

(4) 

cf) 2 e sin cj . !) ^ > 

(5) 
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where 4> x and 0 2 are the eccentric anomalies at times and t 2 , t p is the time 
at pericenter, and 


n = (/x/a 3 ) 1/2 . 


If e t is a unit vector pointing towards periapsis and ? 2 is a unit vector in 
the plane of motion 90° ahead of e 1 in the direction of motion, then for the posi- 
tion vectors r x and r 2 at times t l and t 2 , 

r 2 = a(cos0 2 -e) e 1 + a (l - e 2 ) 1/2 (sin <f> 2 ) 1 2 , 

r 1 = a(cos <p 1 - e) + a(l - e 2 ) 1/2 (sin0 1 )e 2 , 


Substituting these equations in 


r , 


+ r. 


— — f 

2r , 


where the dot indicates scalar product, we have 


c 2 - a 2 (cos <£ 2 - cos 2 + a 2 ( 1 - e 2 ^ (sin 0 2 - sin 


- 4a" 


1 - e 2 cos 2 (^)(^ 1 +0 2 ) sin 2 * 


( 6 ) 


Adding (2) and (3), 


r i + r 2 = 23 


1 - e cos 


(i)K +< £ 2 ) cos (i) • (7) 


Subtracting (4) from (5), 


n { t 2 ~ t i) ~ $2 ~ “ 2e cos (l) (^i +< ^ 2 ) sin ('2')(^2 ~^i) • 
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Equations (6), (7), and (8) determine the three unknowns a, <p 2 - cf> l , and 
e cos (1/2) + (p 2 Y Let 


cos 


(4) ( a + Z 3 ) “ e cos (4 ) (4> x +0 2 ) . 0 S a + /S < 2a , 


( 9 ) 


a - /3 - 4> 2 ~ 4 > 1 ~ 2a m , 0 < a - /3 < 2rr , 


( 10 ) 


where m is the number of complete circuits made by the particle between times 
t j and 1 2 . 

Equations (6), (7), and (8) become 

c/2a = sin^^ (a + /3) sin^ 2 ) ( a “ /3) (H) 

( r l +r 2 )/ 2a = 1 - c Os(4) ( a +/S) c °s (4) O -yfi) (12) 

= lam + a - ft - 2 cos ( 4 ) (a + /3) sin 4 ) (a - /3) . (13) 


The two inequalities in Equations (9) and (10) are geometrically equivalent 
to the shaded region of Figure 1, from which it is evident that 0 < a < 2a and 
- a ^/3 <a. We can also obtain 0 -S a < 2a by adding the inequalities in (9) and 
(10); and if we add /3 - a to each part of inequality (9) and divide the result by 2, 
we obtain - (a - /3)/ 2 i /3 < a + (a - /3)/2 or - a £ /3 < a. 

With appropriate trigonometric identities, Equations (11) and (12) become 


cos /3 - cos a - c/a 


cos /? + cos CL 
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Solving these two equations, 


cos a 1 - s / a 
cos ~ 1 + 2KE 


where we have defined 




Since cos a - 1 - 2 sin 2 (a/2), (14) can be changed to 


E = - sin 2 (a/2) , 0 £. a < 2jt . (16) 

Similarly (15) becomes 

KE = - sin 2 (y8/2) , - 77 < j3 < n (17) 

K = (s-c)/s = ( r i + r 2 -c )/ 2s 

= [(r x + r 2 ) 2 -c 2 ]y4s 2 . 

Introducing (1) we have 



K = 

( r i 

r 2 /2s 2 )(l + 

co s 6 ) 



= 

( r i 

r 2/ /s 2 ] cos 2 

(9/2) . 


Substituting (16) in (17), 





sin (fi/2) 

- q 

sin 

(a/2) , 

- 77 — /3 < 7 T 

(18) 

q - 

± Vic 

= 

[(', n)‘7' 

sj cos (&/2) . 

(19) 

Note that the sign of 

q is taken 

care of by the angle 0: 



1 > q 

> 0 

if 0 < 

• 

$ S n , 



0 > q 

> - 

1 if TT < 

6 < 2tt , 



6 



We can introduce E into (13) since 


n(t 2 -t,) - (,vV) V Mt a -t,) = (-E)«T 


where 


T = (8q/s) 1/2 (t 2 - t t )/ s . 


( 20 ) 


T " (-E) 3/2 211177 + a - f3 - 2 cos 


(- 3 ) (a + /?) sin (-j) (a - /3)J 


( 21 ) 


This can also be written as 


T - (-E) 3/2 [2m7T + a - fi - (sin a - sin/3)] . 


( 22 ) 


Substituting (16) into (22), 


T sin 3 (a/2) 


2mrr + a - fi - sin a + sin/5 . 


(23) 


Equations (18) and (23) with 0 ^ a < 2tt are Lambert's equations for elliptic 
motion. Given T and q, they are to be solved for a and/3, after which it is a 
simple matter to find all other quantities associated with the orbit. 

It is customary in the literature [e.g., 2] to consider T as a function of E (or 
a) and break the elliptic case of Lambert's theorem into four cases, depending 
upon the sign of q and whether from (16) a/ 2 is taken in the first or second 
quadrant. The choice of E as the independent variable makes T a double valued 
function. This problem can be avoided by choosing a as the independent variable. 
However, an even better choice will be discussed in the next section. 
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By a derivation very similar to that for the elliptic case, one finds for the 
hyperbolic case, 


T - - E~ 3/2 [y- 3 - (sinhy- sinh 8)] , (24) 

E = sinh 2 (y/2) , (25) 

sinh (S/2) « qsinh(y/2) . (26) 


When m = 0, Equations (22) and (24) break down for E = 0 and suffer from 
a critical loss of significant digits in the neighborhood of E = 0. To remedy 
this (22) is written in the form 


T = cr(-E) - qKo-(-KE) (27) 

o-(u) = 2 [arcsin u 1/2 - u 1/2 (1 - u) 1/2 ]y/u 3/2 

Replacing arcsin u 1/2 and (l-u) 1/2 by series [3] 
cr(u) = 4/3 


a n = 1*3- 

A similar procedure produces the same series for the hyperbolic case. 

For the parabolic case we have E = 0, in which case the series gives 

T = (4/3) (l~q 3 ) . 

Thus with m = 0 we have a series which is valid for elliptic, hyperbolic, and 
parabolic transfer provided 0< x < 2 1/ 2 . 


00 

n = 1 


U < 1 


5 ••• (2n- 1 )/ 2 n-2 ( 2n + 3 )n ! 
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A UNIFIED FORM OF LAMBERT’S EQUATIONS 


As mentioned in the previous section, T is a single- valued function of a. 
However, a better behaved function is obtained if we choose as the independent 
variable 


x - cos (a/2) , 


-1 < x < 1 , 


- cosh (7/ 2) , x > 1 . 


We then have, for both elliptic and hyperbolic transfer, 

E = x a - 1 . 

For parabolic transfer, it is obvious that we should let x = 1. 
For the elliptic case let 


y - 

sin (a/ 2 ) 


(-E) 

1/2 

z = 

cos (// 2) 

= 

(1 + 

KE) 1/2 

f = 

sin (3) (a 

-/3) 

= 

y(z - qx) 

g = 

cos (3) (a 

~/3) 

= 

xz - qE 

h = 

(l)(sina 

- sin 0 ) 

= y(x-qz) 

k = 

arctan ( f/g) , 


0 — X — v 
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It then follows from (21) for the elliptic case that 

T = 2(m7T + X - h^/y 3 . 

For the hyperbolic case let 


sinh(y/ 2 ) = 

E V2 


cosh (y/ 2 ) - 

(1 + 

KE ) 1/2 

sinh ^7 

*) (7- S) 

- 

y(z - qx) 

cosh ^7 


- 

xz “ qE 


Note that 0 y - 8 < 00 since 0 — f < 00 . Let 

h = [2) (sinhy- sinh 8) = y(x-qz) 


It follows that 


(2)0*- s > 


arctanh ( f/g) 



(i)M<f + g>/(g-f>] 

(i)^ n [c f + g> 2 / (g 2 - f 2 )] 


■in ( f + g) , 
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Thus for the hyperbolic case 


T = 2 [h - ln( f + g)] /y 3 . 


It is now apparent that, given qand x, the following steps produce T for all 
cases: 

1. K = q 2 

2. E = x 2 - 1 

3. P - |e| 

4. If p is near 0, compute Tfrom (27) 

5. y - p 1/2 

6. z = (1 + KE) 1/2 

7. f ■= y(z-qx) 

8. g = xz - qE 

9. If E < 0, A. = arctan ( f/g) ,d = mir + \, 0 < X. < 7r 
IfE> 0,d = ln(f + g) 

10. T = 2(x - qz - d/y )/E. 

The following formula for the derivative holds for all cases except for 
x = 0 with K = 1 and for x = 1. 


dT/dx = (4-4qKx/z- 3xT)/E . 


If xis near 1, differentiate (27) to obtain 


dT/dx = 2x[qK 2 cr' (-KE) - cr' (-E)] , 
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o-' (u) 


n 3 a U 


da / du 



n- X 


The derivative in the ease of x = 0 with k = 1 will be discussed in the next 
section. 


AUXILIARY FORMULAS 

In this section we will show how to obtain a number of useful quantities 
associated with the two-body orbit, assuming Lambert’s problem has been 
solved for x. 

In the derivation of Lambert's equation for the elliptic case, a and /3 are 
defined in such a way that 


4> 2 ~ 4>\ ~ a - /3 + 2mv - 2 (A. +m7T) . 


Using (22) the eccentric anomaly difference can also be written in the form 
<fi 2 - 4> x - (“E) 3/2 T + sin a - sin/3 

= y 3 T + 2y(x- qz) . (28) 

From (28) we have 

sin (4 > 2 - 4> t ) = 2y( z - qx) (xz - qE) 

= 2y(x-qz) +4y 3 q(z~qx) 

4> 2 ~ 4>i ~ sin (4> 2 ~<t> 1 ) = y 3 T - 4y 3 q(z-qx) 

1 - cos ($ 2 - 4> x ) = 2y 2 (z - qx) 2 . 


(29) 

(30) 
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We now obtain a formula for the time derivative r at time t x . Kepler’s 
equation in the elliptic ease can be written in the form [4] 

(^/a 3 ) 3/2 (tj-tJ = 4> 2 ~ 4 > % + r x r x [l - cos (<£ 2 -<£^^<> 3 ) 1/2 


" ( 1_r x/ a ) sin (<£ 2 - <£ x ) . 

Substituting 1/a = 2y 2 /s, t 2 - t x = s 3/2 T/(8/i) 1/2 , and making use of (29) 
and (30) gives 

(2//is) 1/2 (z - qx) r x r x = 2q - (2r x /s) (xz - qE) . 

Multiplying through by z + qx, we have, since ( z - qx) ( z + qx) = 1 - K = c/s 
and (z + qx) (xz - qE) = x + qz, 

(2//is) 1/2 c r x r x = 2 qs(z + qx) - 2r x (x + qz) 

= 2qz(s~r x ) + 2x(Ks~r x ) 

Ks - r x = (l-c/s)s - r x = s - c ~ r x - r 2 - s . 

Thus we have finally that 

r x = (2/is) 1/2 jqz(s - r x ) -x(s - r 2 )Jy/cr x . (31) 

At time 1 2 we find 

r 2 = ( 2/xs ) 1 7 2 jx (s - r x ) - qz (s - r 2 )J /cr 2 , (32) 
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By a similar procedure we can show that (31) and (32) hold also for the 
hyperbolic and parabolic cases. 

Having x we can find the semimajor axis a or its reciprocal 

1/a = 2y 2 /s . 


We know that 

e cos <f> = 1 - r/a , 

esinc/) = rr/(/i.a) 1/2 , 

where r is the magnitude of the position vector at time t. 

Thus for the eccentricity we have 

e 2 = ( 1 - r/a) 2 + (rr ) 2 //xa . 

For the semilatus rectum p we have 

p = a(l-e 2 ] = 2r - r 2 /a ~ (rr) 2 //x 

and for the value of r at the point of closest approach to the center of attraction 
we have 

P 

r p - 1 + e 


For the speed we have 



For the component of velocity perpendicular to a radius vector, 


v. = <>p) 1/2 /r . 


The above formulas hold for any type of two-body motion. Either of the 
subscripts 1 or 2 can be placed on r, r, and v. 

If r j < 0 and r 2 > 0 or if fj and r 2 have the same sign with 6> n , pe'riapsis 
passage will occur between times t x and t 2 , in which case it may be of interest 
to compute r p . 

If the plane of motion is known, the velocity v at either time 1 1 or time 1 2 
can be written in any convenient coordinate system, since the components r and 
v 9 are known. The plane of motion can be found from the position vectors r, 
and r 2 at times 1 1 and t 2 , provided they are not parallel. Since this is a com- 
mon case, we will express v 1 and v 2 (the velocities at t x and t 2 ) in terms of 
r j and r 2 . 

For the velocity v t we have 


V r 1 + 


'61 


where v rl is along r x and v d x is in the plane of motion perpendicular to r j and 
in the direction of motion, i.e., in the direction of increasing true anomaly. We 
have 


? rl " (V r l)' 7 l’ 

Vei = C l^l +C 2 7 


where Cj and c 2 are to be determined. 


’ %i * 0 " c i r i + C 2 7 l - 7 2 

r 2 * v 01 = r 2 Vq j sin 9 = c 1 ? 1 • r 2 + c 2 r 2 2 . 


15 



Since r x • r 2 - r x r 2 cos 8 we have 


r i c i + ( r 2 cos 6 ) c 2 " 0 

^r 1 cos0jc 1 + r 2 c 2 - v^sin# , 

Solving for Cj and c 2 we have 

" {^i- v ei cotd )Pi/ T i) + (v gl csc0)(r 2 /r 2 ) . 

In a similar way we find 

^2 = ~( v 02 csce )Pi/ r l) + C r 2 +v e2 cot0 )(*2/ r 2) ■ 


Figure 2 shows T as a function of x. Note the discontinuity in the slope for 
x = 0 with K - 1, For K = 1, z = | x|. Thus we are led to consider four cases: 
q = ±1 with x>0, q = ±1 with x < 0. Examination of the formulas for dT/dx 
in these cases reveals that if q = 1 we have a left-hand derivative of -8 and a 
right-hand derivative of 0 at x = 0. If q = 1 we have a left-hand derivative of 
0 and a right-hand derivative of -8 at x = 0, 


With further analysis we find that the cases (q = 1, x <0) and (q = -1, x >0) 
represent rectilinear orbits. 

For m = 0, T is a monotone function of x, making possible a simple numeri- 
cal procedure for solving Lambert's problem. Figure 2 is for the elliptic case, 
Figure 3 for the hyperbolic case, the parabolic ease occurring at x = 1 in both 
figures. Figure 4 shows a small region of Figure 3 where d 2 T/dx 2 is negative. 
If the Newton-Raphson method is being used to find x, a switch should be made 
in this region to the secant (regula falsi) method. 

No solutions of Lambert's problem exist in the shaded regions of Figures 2 
and 3. x = 1 (m>0) and x = -1 are vertical asymptotes. T-o asx-®. 
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